Renormalization effects in interacting quantum dots coupled to superconducting leads 
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We study subgap transport through an interacting quantum dot tunnel coupled to one normal 
and two superconducting leads. To check the reliability of an approximation of an infinitely- large 
gap A in the superconducting leads and weak tunnel coupling to the normal lead, we perform a 1/A 
expansion, and we analyze next-to-leading order corrections in the tunnel coupling to the normal 
lead. Furthermore, we propose a resummation approach to calculate the Andreev bound states for 
finite A. The results are substantially more accurate than those obtained by mean- field treatments 
and favorably compare with numerical exact results. 



PACS numbers: 74.45. +c,73.63.Kv 

I. INTRODUCTION 

The tunnel coupling to superconducting reservoirs can 
induce superconducting correlations in quantum dotsP^ 
Interesting linear- and nonlinear transport effects aris- 
ing from Andreev reflection 3 -^ can be studied in hy- 
brid normal-superconducting systems comprising quan- 
tum dots in the presence of strong Coulomb interac- 
tion. During the past years, there has been tremen- 
dous experimental progress in coupling quantum dots 
to superconductors pHSU p rom i nen t examples include the 
observation of — tt transitions in a Josephson junction^ 
and splitting of Cooper pairsP^H^l 

An exact theoretical treatment of the transport prop- 
erties of such devices requires to account for Coulomb 
interaction, superconducting correlations, and non- 
equilibrium at the same time. To circumvent this chal- 
lenge, various approximation schemes have been pro- 
posed. In the case of vanishing Coulomb interaction, 
an exact solution can be obtained within a scattering 
approach. 23 As a first step beyond, Coulomb interaction 
has been treated perturbativelyf^HHI which, however, is 
not justified for the large charging energies that are typ- 
ical for small quantum dots. An alternative possibility 
is to allow for strong Coulomb interaction but perform 
a perturbation expansion in the strength of the tunnel 
coupling between dot and leadsPSEB However, in the 
limit of an infinitely-large superconducting gap, an exact 
treatment of both the Coulomb interaction and the tun- 
nel coupling between quantum dot and superconducting 
leads is possible.^ES! i n this case, quasiparticle tunnel- 
ing is completely suppressed and transport from and to 
the superconducting leads is fully sustained by Cooper 
pairs. The formation of Andreev-bound states (ABSs) 
in the quantum dot indicates that superconducting cor- 
relations are induced via the proximity to superconduct- 
ing leads. The dependence of the ABS energies on the 
quantum dot's level position is directly reflected in the 
transport spectrum of the system. 

A diagrammatic real-time approach to transport 
through quantum dots coupled both to normal and su- 



perconducting leads has been introduced in Ref. |3T| This 
method has been applied in the limit of an infinitely- 
large superconducting gap, A — >> oo, to study a va- 
riety of transport phenomena in different quantum- 
dot setups) 36 " 43 * including, for example, pure spin- 
current gene ration ,^ shot-noise suppression,^ Cooper- 
pair splitting) 39140 * and time dependent driving.^1 

The calculations in the infinite-gap limit neglect quasi- 
particle contributions to transport. They, furthermore, 
approximate the ABS energies since the latter are a func- 
tion of the superconducting gap. In the experiments, 
the superconducting gap A is often of the same order of 
magnitude as other energy scales such as charging energy 
or tunnel-coupling strength. For this reason, it is quite 
natural to investigate the quality of the approximation 
introduced by considering the infinite-gap limit. This is 
the main goal of the present paper. Furthermore, we an- 
alyze the effect of next-to-lowest order corrections in the 
tunnel coupling to the normal leads. 

A finite superconducting gap in the leads affects the 
Andreev and Josephson transport through quantum dots 
in two ways. First, it modifies the ABSs. This affects the 
position of bias or gate voltages at which a new transport 
channel is opened. Second, the value of the current away 
from these threshold voltages is changed. To address the 
effect of a finite A, we start by performing a systematic 
expansion in 1/A around A —> oo. Exemplarily, we will 
focus on a quantum dot tunnel coupled to one normal and 
two superconducting reservoirs and compare the currents 
obtained by the 1/A expansion with the A — » oo limit 
We find that away from the threshold voltages, both the 
Andreev and the Josephson currents are hardly modified, 
i.e., the A — » oo limit provides quite an accurate descrip- 
tion of subgap transport. The 1/A expansion, further- 
more, indicates a renormalization of the ABS energies. It 
can, however, only predict the direction towards which 
the energies are shifted. To determine the position of 
the ABS, we go beyond the 1/A expansion and perform 
a partial resummation of diagrams. The results within 
this approach are considerably more accurate than those 
obtained by a simple Hartree-Fock treatment and com- 
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FIG. 1: (Color online) Setup: a quantum dot is tunnel coupled 
to one normal and two superconducting leads. 



pare favorably with the numerical renormalization group 
(NRG) data of Martm-Rodero and Levy Yeyati!^ 

Finally, we present (for the A — > oo limit) results from 
a systematic perturbation expansion to next-to-lowest 
order in the tunnel coupling to the normal leads. In 
this way, we go beyond the regime of we ak tunnel co u- 
pling but do not cover Kondo correlations! 15 ' 16 1 24 ^ 27 1 45 ^ 
We find that the next-to-leading correction leads to a 
renormalization of the gate-voltage position at which a 
— 7r transition occurs in the Josephson current. This 
renormalization also affects the Andreev current and the 
average quantum-dot charge. 



II. MODEL AND METHOD 

A. Hamiltonian 

We focus on a three-terminal setup composed of a nor- 
mal conductor and two superconductors tunnel coupled 
to an interacting single- level quantum dot, see Fig. [I] 
The system's Hamiltonian is given by H — H^ ot + 
^Z^Hrj + H tunn ^, with r] G {L,R, N}. We assume the 
normal lead to be a reservoir of noninteracting electrons 
and model the superconductors with the mean- field BCS 
Hamiltonian 



where e is the energy of the spin-degenerate single- 
particle level, and U the on-site Coulomb repulsion. The 
dot's annihilation (creation) operators of spin <j are given 

by da \ and n a = d\d a is the corresponding number op- 
erator. Tunneling between dot and leads is described by 
the tunneling Hamiltonian 



tftunn,^ =Vr ] J2 ( C lk<r d * + H ' C ') 
ka 



(3) 



with the spin and momentum independent tunnel matrix 
elements V^. We define the tunnel-coupling strengths as 
= 27r\V r] \ 2 p r] , where p v is the density of states in lead 
rj that we assume to be independent of spin and energy. 

B. Diagrammatic Technique 

1. Effective Dot Hamiltonian 



The main idea of the diagrammatic real-time technique 
is to integrate out the leads' degrees of freedom to arrive 
at a reduced density matrix for the quantum dot. As 
a basis for the reduced density matrix one may use the 
eigenstates of the decoupled quantum dot, iifdot- This 
is a convenient choice for weak tunnel couplings to the 
superconductors!^ In the opposite limit of strong cou- 
pling to the superconductor, a resummation of an infi- 
nite number of diagrams is required. This resummation 
can be carried out exactly in the A — >• oo limit M*^ For 
a systematic 1/A-expansion, however, it is more conve- 
nient to work in the eigenbasis of H e fi = ^dot — H p , 
with the pairing Hamiltonian H p = x*d^d^ + x^j^j an d 

tonian 



jt]=l r r^e^ 77 . The remaining part of the Hamil- 
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(Ae^ 77 c^-k^c^ + H.c.) , i s treated diagrammatically as a perturbation. Thereby, 



(i) 

with e^k being the single-particle energies and c^ a are 
the annihilation (creation) operators for electrons in lead 
r] with momentum k and spin a. Here, A is the modulus 
of the superconducting pair potential, which we assume 
to be the same in both superconductors, and ^ is the 
corresponding phase in lead rj = S with S G {L, R}. We 
measure all occurring energies with respect to the Fermi 
level of the superconductors, i.e. ^l,r = 0. Furthermore, 
we assume all occurring excitation energies of the dot to 
lie inside the superconducting gap of the leads so that 
only subgap transport takes place. 

The Anderson Hamiltonian describing the single-level 
quantum dot reads 



H p has been chosen such that in the A — >• oo limit the 
diagrammatic contributions stemming from H p exactly 
cancel those from the tunneling to the superconducting 
leads, i.e., H e ^ describes the hybrid system of quantum 
dot and superconducting leads in the limit of an infinite 
superconducting gap. 



The 



eigenstates of H e R 
Te-^/yiT^lO)- 
eigenenergies e 
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ire |t), 



ID, 



1 ± 



2e A 



and 

~\d) 



with 



e, and E± = | ± e^, respectively. Here, 



#dot = y^edldg- + Un t n^ 



(2) 



<D X = arg(x), e A = ^(f ) 2 + | X | 2 , and S = 2e + U is 
the detuning between the energies for empty and doubly- 
occupied quantum dot. The superconducting proximity 
effect is indicated by the fact that |0) and \d) are no 
longer eigenstates but appear as linear combinations in 
|d=). The mixing between |0) and \d) becomes largest 
around zero detuning, S ~ 0. 
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The excitation energies associated with H e ^ are given 
by the differences of the eigenenergies of states with even 
and odd dot occupation numbers, 



rpA— >oo 



= 7 + 7^A, 



(5) 



with 7', 7 = ±. They are nothing but the Andreev bound 
state (ABS) energies^ in the A — >• oo limit. 



2. Generalized Master Equation 

The system's dynamics is determined by a general- 
ized master equation for the reduced density matrix p rec b 
which is obtained by integrating out the lead's degrees of 
freedom. Its matrix elements Pg* = (X1I/WIX2) obey, 
in the stationary limit, the generalized master equation 



i(E x 



E ) P Xl 



XiX2 



X2X2 X 2 ' 



(6) 



with W XXl Ir] being the generalized current rates. For a 

XX2 

systematic perturbation expansion in powers of Tn and 
1/A one needs to expand P,W, and W Ir] in Eqs. Q and 

©• 

We refer to the current flowing out of the normal 
lead as the Andreev current and to the difference of the 
currents between left and right superconductor as the 
Josephson current. In the limit A — » 00, the Andreev 
current and the Josephson current are directly connected 
to the isospin visP^ 
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(10) 

(ii) 



For the 1/A corrections to the currents, however, one 
needs to use the more general Eq. (p)]). 



4- Diagrammatic Rules 



where W x ^ are generalized rates that can be computed 



in a diagrammatic way, see Sec. II B 4 As explained in 
the previous section, we work in the basis {t?b+? — }■ 
This has the consequence that the off-diagonal matrix 
elements on the right hand side of Eq. ^ vanish to lowest 
order in T^r and 1/A as long as ca ^T^. Only their 
next-order correction remains finite. 

Superconducting correlations in the quantum dot are 
associated with coherent superpositions of an empty and 
a doubly-occupied dot. This motivates the following def- 
inition of an isospinj 30 ! 31 ! 



p d + p° 



r>d r>0 

; ~ r d 



Po 



(7) 



formulated in the basis {t, i, 0, d} (the transformation to 
the basis {t, i, +, — } is straightforward). A finite value of 
the isospin components I x and/or I y indicates the pres- 
ence of superconducting proximity effect. 

It turns out that the generalized master equations for 
the isospin components do not couple to the probabilities 
of the quantum dot to be singly occupied with either spin. 
They can be written in the form of a Bloch-like equation 







dl 

dt 



A-RI + IxB, 



(8) 



where the terms with A, R, and B describe the genera- 
tion, relaxation, and coherent rotation of the isospin. 



3. Current formulae 



The current in lead 77 is given by 



V w xx , 

Z-^ xx 2 



pXi 

x' 2 ' 



(9) 



Starting from Refs. [31] and [37] we derive the rules for 
calculating the generalized rates and the generalized 
current rates in the basis of the effective Hamiltonian. 
The diagrammatic expansion is performed with respect 
to the tunneling parts of the Hamiltonian and, in 
addition, Hp. In the following, we assume the two 
superconductors being at zero chemical potential, i.e. 
Ml = I^r = and a symmetric tunneling coupling, 

r L = r R = r s . 

(1) Draw all topologically different diagrams with 
fixed ordering of the vertices with respect to (real) time. 
The vertices may be bullets (representing tunneling) or 
crosses (for H p ). The bullets are connected in pairs by 
tunneling lines carrying energy uji. The tunneling lines 
can be normal or anomalous. For each anomalous line, 
choose the direction (forward or backward with respect 
to the Keldysh contour) arbitrarily. 

(2) For each vertical cut between two vertices, assign a 
factor l/(A£? + iO + ), where AE is the difference between 
the left-going and the right-going energies, including the 
energies of the dot states E x and the tunneling lines uj{. 

(3) For each tunneling line arising from the normal lead 
assign a factor ^^n/n^n) and for each tunneling 
line arising from a superconducting lead assign a factor 



1^1 



r 6(K|-A)/±(a;, 



where f£(u) 



XX1X2 



27r ^vWp- a2 

f + (uj-fi N ) = f(uj-ii N ) = {1 + exp [(u - fj, N ) /ksT}}- 1 
and fx(u) = f~(u- Mat) = 1 - f(u - Mat). The upper 
(lower) sign applies for lines going backward (forward) 
with respect to the Keldysh contour. For anomalous 
lines, multiply an additional factor ±sign (c^) 

Moreover, assign a factor e l ® x e - *^ 77 for an outgoing 
and e~ z ® x e 2$r? for an incoming anomalous line, where 
$ x = arg ( X ) and x = h r s E ne s ^ • 
(4) Each bullet vertex connects a dot state |a), 



with a G {+,—}, to another dot state \a), with 
a G {t?l}> an d to a tunneling line describing a 
tunneling electron of spin a' G {t?^}- For each 

bullet vertex assign a factor ^l-S a S a S a ,^/V2, 

where S a= + = S a =^ = S a >=-f = 1 and 

S a =- = Sa=l = Sa'=l = ~ 1- 

(5) Each cross vertex connects a state \a) with 
a G {+,— } to a state \a!) with a' G {+,—}• 
For each cross vertex with a = a' assign a factor 



(-l)SaJl 



c5 2 



|x| and for a ^ a' assign a factor 



6 \xV 



2e A 

(6) Assign an overall prefactor —i. 
Furthermore, assign a factor —1 for each 

(a) vertex on the lower propagator; 

(b) crossing of tunneling lines; 

(c) outgoing (incoming) bullet vertex ending (starting) 
in the state |+) or | "f); 

(d) outgoing (incoming) anomalous tunneling line, in 
which the earlier (later) tunnel vertex with respect to 
the Keldysh contour is connected to a tunneling spin-up 
electron. 

(7) For each tunneling line, integrate over the energy Ui. 
Sum over all diagrams. 

The generalized current rates W xx ?~ Ir] are evalu- 
ated in the following way: 

(8) Multiply each diagram that contributes to the 

corresponding generalized rate W**, 1 and where the 

XX2 

rightmost line is associated with lead 77 with a factor: 

(a) for the rightmost line being a normal line: 1 if the 
line is going from the lower to the upper propagator, -1 
if it is going from the upper to the lower propagator, 
and otherwise; 

(b) for the rightmost line being an anomalous line: 1 
for incoming lines within the upper and outgoing lines 
within the lower propagator, —1 for outgoing lines 
within the upper and incoming lines within the lower 
propagator, and otherwise. 



We want to remark that all terms arising from rule 
(5) are exactly cancelled by the terms arising from the 
superconducting tunneling lines in the infinite-gap limit. 
This reflects the fact that the effective Hamiltonian com- 
pletely describes the hybrid system of quantum dot and 
superconducting leads in the infinite-gap limit. To illus- 
trate how the graphical representation is translated into 
analytical expressions we show two exemplary diagrams 
here. The first example is a diagram containing one line 
associated with the normal lead and one line associated 
with the left superconductor, see Fig. [2] (a), and it reads 




FIG. 2: Example diagrams contributing to the rate W^. 



4 1^1 ^Ws^-'^tiMf-M 



1 



uj n + E- - e + i0+ ujn + ojl + ^0+ ujn + E+ — e + i0+ 

The second example is a diagram, see Fig.[2](b), contain- 
ing a line associated with the right superconductor and 
a cross vertex and it yields the following expression 



1 I ^Tse^-^f-^n) 



2tt 



E + -ojr- e + i0+ 

s\x\ 



sign(w fl ) 



E- - lor - e + i0+ 
A 



x e (M -A)^Ji-^. 



III. FINITE-GAP EFFECTS 
A. 1/A Expansion 

We assume symmetric tunnel couplings to the super- 
conducting leads, = Tr = Ts- Furthermore, we 
choose a gauge such that = — &r = 3>/2. For both 
the Josephson and the Andreev current we calculate the 
value for the A — >• 00 limit and the correction to first 
order in 1/A. The tunnel coupling to the normal lead, 
Tat is accounted for in lowest non- vanishing order, i.e., to 
zeroth order for the Josephson current and to first order 
for the Andreev current. 

Figures [3] and [4] show the Josephson current and the 
Andreev current in the limit of A — 00 as a function 
of the level position e and the chemical potential of the 
normal lead [in. This is the limit considered in Ref. [31] 
and is used as a reference for the present discussion. A 
cut through the density plots at a fixed level position e = 
—0.4/7 is shown in Fig. [5] (solid blue line). The currents 
display a step-like behavior with the position of the steps 
reflecting the ABS energies. This is compared with the 
currents expanded up to first order in 1/A, dashed red 
lines in Fig. [5] From this comparison, we can draw two 
conclusions. 

First, away from the ABS, the solid and dashed curves 
almost coincide, i.e., the 1/A correction to the current 
is small. This leads to the conclusion that the A — >• 00 
calculations provide a useful and quantitatively accurate 
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FIG. 3: (Color online) Density plot of the Josephson current 
in the limit of A — >• oo as a function of the level position e 
and the chemical potential of the normal lead fiN- The other 
parameters are F s = 0.2U, F N = 0.001?7, = 7r/2, and 
k B T = 0MU. 



FIG. 4: (Color online) Density plot of the Andreev current 
in the limit of A oo as a function of the level position e 
and the chemical potential of the normal lead fiN- The other 
parameters are F s = 0.2U, F N = 0.001?7, = 7r/2, and 
k B T = 0.01U. 



description of sub-gap transport also for finite values of 
the superconducting gap. 

Second, we find peaks at the positions of the ABSs in 
the 1/A correction to the current. While the peaks them- 
selves are artefacts of the 1/A expansion, they carry im- 
portant physical information about the renormalization 
of the ABSs as compared to their A — >• oo value. This 
can be understood in the following way. The currents 
vary stepwise at the position of the ABSs. The latter are 
renormalized due to finite A. Therefore, a systematic 
expansion in the small parameter with which the renor- 
malization scales gives rise to the occurrence of peaks at 
the original position of the steps. The sign of the peaks 
relative to the sign of the steps indicates the direction of 
the renormalization. The ABS is renormalized towards 
(away from) zero if step and peak have the same (the op- 
posite) sign. In the examples plotted in Fig. [5] the curves 
indicate a renormalization of all four ABSs towards zero. 



B. Renormalization of Andreev Bound States 

The 1/A expansion indicates that the ABSs are renor- 
malized as compared to the A — >• oo limit. It even pre- 
dicts the sign of the renormalization. It does, however, 
not provide a possibility to calculate the ABSs at finite 
A. For this task, alternative approaches are required. In 
the limit of vanishing Coulomb interaction, the problem 
can be solved exact lypHHl A mean- field treatment of the 
Coulo mb in teraction provides easily an estimate of the 
ABSsP233 Numerical renormalization group (NRG) cal- 



culations are exact but numerically demanding! 44 * 58 ! ^ 
Instead, we propose a resummation approach that is 
more accurate than mean field but less computationally 
heavy than NRG. 

In the following, we consider the limit of weak tunnel 
coupling to the normal lead, i.e., we discuss the spectrum 
to zeroth order in Fjy. Furthermore, we choose = 
without loss of generality (to consider the case of a finite 

one has to replace Fs by Fs cos(<£/2); to model the 
case of a single superconducting lead, Fs needs to be 
replaced by F$/2)- 



1. Exact Green's function for the noninteracting dot 

For a noninteracting dot, U = 0, the exact Green's 
function for spin a is given, in Nambu space, by 

G r = (oj • I 2x2 - h a (U = 0) - Es)" 1 , (12) 

with I nxn being the identity matrix in n dimensions. The 
matrix 



K(U = 0) 



e CT 
-e_ 



(13) 



accounts for the single-particle energies, and the self- 
energy due to tunneling reads 



(14) 
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FIG. 5: (Color online) Plot of (a) the Josephson current and 
(b) the Andreev current as a function the chemical potential 
of the normal lead /xjv- The solid (blue) curve shows the 
A — » oo limit, the dashed (red) curve includes the first l/A 
correction. The other parameters are e = — 0AU, Fs = 0.2C7, 
T N = 0.001C7, A = W , $ = tt/2, and fc B T = 0.0117. 
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FIG. 6: (Color online) Plot of the U = Andreev-bound 
states for a tunnel-coupling strength to the superconducting 
lead of Fs = 0.25A as a function of the level position e. All 
energies are normalized to A. The solid (blue) line shows 
the exact result, the dashed (black) line the bound states 
computed by means of the A — » 00 expression, which in the 
non-interacting case is simply ±6a- 



The ABSs, i.e. the excitation energies, are found as the 
real part of the poles of the Green's function's determi- 
nant. An example is shown in Fig. [6j In the noninter- 
acting case, there are only two ABSs. Their energies for 
finite A are reduced as compared to the A — » 00 limit. 
In particular, the curvature of the energy as a function of 
the bare level position e changes for large |e|, such that 
the ABSs lie within the window [—A, A]. 



2. Hartree-Fock approximation 

A Hartree-Fock (HF) treatment of the Coulomb in- 
teraction has been used in Ref. |2l] to fit the experi- 
mental data. This HF approach corresponds to replace 
h a (U = 0) in Eq. O by 



Kill) 









-e_ 



(15) 



where eh 
Eq. @ 



e^ + L7. For U = 0, the exact Green's function 
is recovered. Furthermore, it reproduces the 
exact ABSs in the A — >> 00 limit. 

An improved version is self-consistent HF, in which 
the dot occupations are calculated self-consistently.^ In 
the following, however, we use the simpler version of HF 
for comparison with the results from the resummation 
approach. 



3. Resummation Approach 

An important feature of the diagrammatic technique 
presented in this paper is that Coulomb interaction can 
be included beyond the mean-field level. The downside 
is that for an exact solution one has to sum up infinitely 
many diagrams. A great simplification is achieved by 
working in the basis {t, |, +, — }. As discussed above, 
in the A —> 00 limit, all diagrams that contain super- 
conducting lines are cancelled by cross vertices. This 
means that all contributing diagrams contain normal tun- 
nel lines only. 

At finite A, the situation is different. Now, both su- 
perconducting lines and cross vertices dress the original 
normal tunnel lines. We define the propagator FI(uj) as 
the diagram part between two vertices that are connected 
by a normal tunnel line with energy uj running from right 
to left. This propagator obeys a Dyson equation 



(16) 



where U^°\uj) is the free propagator (without supercon- 
ducting lines and cross vertices) and E(co>) the self energy. 

In the A — >• 00 limit, the self energy is exactly zero. 
For an approximative treatment of the case of finite A, 
we now include all self energy parts that contain one su- 
perconducting tunnel line or one cross vertex. Examples 
are shown in Fig. [7[ 

Since the Hilbert space is four dimensional, the di- 
mension of the corresponding Liouville space and, thus, 
of the matrices appearing in Eq. (16) is, in general, 16. 



However, as a consequence of spin and charge conserva- 
tion while tunneling, the Dyson equation decouples into 
blocks of 4 x 4 matrices. The presence of the normal 
tunneling line running from right to left which carries a 
definite charge and spin puts constraints on the possible 
states at the beginning and end of the upper (xi) and 
lower (X2) contour. In the following, we choose, without 
loss of generality, (xi,x 2 ) e {(-, t), (i, +), U, -), (+, t 
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(a) 



(6) 





FIG. 7: Examples for diagrams that contribute to the self- 
energy. 



)}, which is the relevant case when the normal tunneling 
line carries spin ^. Then, 

(n(°)H) _1 = 

uj ■ I 4x4 - dia,g(E A ,+,-, E A ,-,-, E A ,-,+ , E A ,+,+§17) 
for the inverse free propagator, and the self energy 

(W^~{uj) W~^(oj) W~j(uj) W^ + (u)\ 
W$(w) Wl l + (u) wiUcv) W#(u) 
W^~(uj) Wi l + (uj) Wi l _(oj) W^(uj) 



(18) 

is given by the generalized rates W^lxl' ( u ) tna ^ contain 
a left-going external line with energy uj. 



We solve Eq. (16) for the full propagator to obtain 

-l 



(n<°>M) '-em 



(19) 



The excitation energies of the proximized dot are probed 
by the external energy uj and they are given by the real 
poles of the full propagator. As mentioned above, we 
approximate the self-energy E(cj) by including only di- 
agrams that contain a single superconducting tunneling 
line or one cross vertex, in addition to the external line 
of energy uj, see Fig. [7| 

Although the resummation scheme does not define a 
controlled approximation, its results are with a few lim- 
itations quantitatively accurate, as we discuss in the fol- 
lowing. It is clear that, by construction, the A — >> oo limit 
is reproduced exactly. Also for the noninteracting case, 
[7 = 0, the resummation approach yields the exact solu- 
tion. As an artifact, however, the resummation scheme 
produces for U = two extra, unphysical, solutions in 
addition to the two correct ones shown in Fig. [6j 

The most interesting regime, however, is the case of 
finite U and finite A. Results obtained from the resum- 
mation approach for different values of T$ and U are pre- 
sented in Fig. [8] as solid (blue) curves. Depending on the 
values of the parameters of the model, we find either two 




FIG. 8: (Color online) Plot of the Andreev-bound states from 
the resummation approach, NRG 62 , and the Hartree-Fock ap- 
proximation as a function of the level position e for (a) U = A 
and T s = 0.5A, (b) U = A and T s = 0.2A, and (c) U = 4A 
and Fs = 0.2A. 



or four solutions. The ABSs are renormalized as com- 
pared to the A — >• oo solution. The sign of the renomal- 
ization agrees with the prediction of the 1/ A expansion of 
Sec. Ill A For comparison, we also show the results from 
HF as dashed (black) lines. While qualitatively similar, 
there are substantial quantitative deviations. First, HF 
seems to underestimate the renormalization of the ABSs. 
Second, the position of the crossing points of the ABS en- 
ergies in the resummation scheme is shifted as compared 
to the HF and the A — » oo result, see Fig. |8jc). 

A reliable quality check of the proposed approximation, 
however, is only given by the comparison with the full 



8 



NRG results.^ We find a remarkably good agreement be- 
tween NRG data, dotted (red) lines in Fig. |8| and resum- 
mation approach. The ABS spectrum is quantitatively 
reproduced in almost the entire parameter space. The 
main difference is that the resummation scheme some- 
times yields four solutions when NRG only predicts two. 
Furthermore, sometimes the resummation approach pro- 
duces solutions with finite imaginary part. Since the 
ABSs in the absence of coupling to a normal reservoir are 
sharp due to the gap in the quasiparticle spectrum, these 
solutions are discarded. This gives rise to the small gaps 
in the solid (blue) line in Fig.|8ja) and (b). Both features 
(extra pair of solutions and poles with finite imaginary 
part) are artifacts of the approximation employed. 



IV. BEYOND WEAK COUPLING TO THE 
NORMAL LEAD 

In the previous section we have discussed how a finite 
superconducting gap in the leads renormalizes the ABSs 
of a quantum dot weakly tunnel-coupled to a normal 
lead. Now, we want to address the influence of the tun- 
nel coupling to the normal lead beyond the weak-coupling 
limit. For this, we go back to A — > oo but include next- 
to- leading order corrections in T^. To distinguish the 
contributions of different order in T^, we introduce in 
the following an index (n) for n-th order. We study the 
Josephson and the Andreev current as well as the aver- 
age charge of the quantum dot. Within a perturbation 
expansion, they start to zeroth, first, and zeroth order, 
respectively, i.e., the next-to- leading order corrections are 
of first, second, and first order, respectively. To evaluate 
them we have to include diagrams with two tunneling 
lines from the normal lead. This increases the number 
of diagrams considerably and we generate and evaluate 
them by means of a computer code. 

It is natural for a systematic perturbation expansion 
in a small parameter that higher-order corrections are 
quantitatively small. Therefore, we focus, in the follow- 
ing, on qualitative features that appear in the region close 
to the point of zero detuning S — 2e + U = when going 
beyond the weak-coupling limit. Around zero detuning, 
the weights of the states |0) and \d) entering the coher- 
ent superpositions |d=) become equal, indicating a large 
proximization of the quantum dot. For bias voltages large 
enough such that the quantum dot has a finite probabil- 
ity to be either empty or doubly occupied, Josephson and 
Andreev currents set in. 

In the regime of large bias voltage, the Josephson cur- 
rent changes sign as function of the quantum dot's level 
position e. This indicates a transition from a 0- to a ir- 
junction behavior. In the limit A oo and weak tunnel 
coupling to the normal leadj^ this — tt transition hap- 
pens exactly at zero detuning, e = — U/2, and then bends 
when the bias voltage approaches the highest ABS, see 
Fig. [3] When including the next-order correction in T^r, 
the position of this — ir transition becomes renormal- 
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FIG. 9: (Color online) Plot of (a) the position of the — 
7r transition of the Josephson current and (b) the extreme 
values of the Andreev current. The chosen parameters are 
k B T/U = 0.05, Y N /U = 0.05, Y s /U = 0.2, and $ = tt/2. 
The thin black lines show the position of the ABS. 



ized, Fig. [9] (a). The origin of this renormalization can be 
understood in terms of an effective, tunnel-coupling in- 
duced field that acts on the isospin introduced in Eq. ^ . 
A finite bias voltage leads to the generation of the z- 
component of the isospin. The tunnel coupling to the 
superconductors leads to an effective field along the in- 
direction. This causes the isospin to rotate and acquire a 
finite ^-component. According to Eq. (10), however, a fi- 



nite x-component of the isospin is needed for a Josephson 
current. This is accomplished with an additional finite 
z-component of the effective field. To lowest order, this 
component is given by the detuning, B^ = J, which ex- 
plains the position of the — ir transition in Fig. [3] The 
next-order correction, renormalizes this position. 

We remark that this renormalization of the — tt tran- 
sition has already been predicted in Ref. [30j in which the 
tunnel couplings to both the superconductors and the 
normal leads as well as the detuning S have been simulta- 
neously treated as small parameters, i.e., the systematic 
perturbation expansion included already both B^ and 
B^ '. On the other hand, the formation of ABSs could 
not be treated in that expansion. In Ref.|3TJ on the other 
hand, no constraints were put on Ts and 5, i.e., the ABSs 
could be described, but the renormalization B^ did not 
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FIG. 10: (Color online) Plot of the average dot charge. The 
chosen parameters are ksT/U = 0.05, Tn/U = 0.05, Ys/U = 
0.2, and $ = tt/2. 



contribute to lowest order in Tat. Only the present cal- 
culation, with arbitrary Ts and S and next-to- leading or- 
der on Tat enables us to address both the ABSs and the 
renormalization of the position of the — tt transition. 

While the Josephson current changes sign at small de- 
tuning, the Andreev current becomes extremal as a func- 
tion of e, see Fig. [4] In Fig. [9] (b) we show the position 
of this extremum in the e-/iiv-plane. (We do not show 
the position of the extremum in the Coulomb-blockade 
regime). To lowest order and for large bias voltage, the 
extremum is at zero detuning. The next-order correction, 
however, leads to a renormalization of this position, for 
the same reason as the renormalization of the — tt tran- 
sition in the Josephson current. 

Finally, we address the average quantum dot charge 
Q, which is related to the z-component of the isospin 
via Q = — e (1 + 21 z ). Again, we focus on the region of 
small detuning where the proximity effect is most pro- 
nounced. In the Coulomb-blockade regime (at small bias 
voltage), the dot is preferably singly occupied, Q — — e. 
A larger bias voltage has the tendency to inject a sec- 
ond electron into the quantum dot. Superconducting 
correlations, however, mix the states of empty and dou- 
ble occupation and, thus, reduce the average number of 
electrons. For zero detuning 5 and to zeroth order in 
Tat, the charge Q remains — e even for large bias voltage. 
Only when increasing \S\ the number of charges increases 
above 1. In the intermediate-bias regime for finite de- 
tuning values the average charge tends towards or — 2e, 
depending on the sign of the detuning. This effect can 
easily be understood by looking at the eigenstates |±): 
in the intermediate-bias regime for positive values of /i/v, 
i.e. Ea,+,- < Mat < £U,+,+ , the state |— ) can already 
be excited while state |+) is still inaccessible. From the 
definition of state | — ) it follows that a positive (nega- 
tive) detuning causes the weight of the zero component 
to dominate over (be dominated by) the double compo- 
nent resulting in an decrease (increase) of the average 
number of electrons. 

1, 



or equivalently I z = 0, for both zeroth and first or- 
der. (We do not plot the position in the Coulomb- 
blockade regime.) As discussed above, the position of 
Q(°V( — e ) — 1 is at zero detuning in the large-bias 
regime and splits in the intermediate-bias regime, with 
Q (0) /(-e) < 1 between and Q (0) /(-e) > 1 outside the 
two lines. When including the next-order correction, the 
situation changes. Now, the line of [Q^+Q^]/(— e) = 1 
is split even for the large-bias regime, opening a region 
with [Q(°) +Q( 1 )]/(-e) < 1 in between. To understand 
this, we analyze I z at zero detuning, 6 = 0. Solving 
the generalized master equatior i 31 * 63 ! immediately yields 
I { z ] (6 = 0) = and I z 1] (6 = 0) = -A^/B^\ where 
=2\ X \ and 




S 2 



7,7 / = ± 



1 2 27rk B T 



In Fig. 10 we show the position of Q/(—e) 



Here W XlX ^ are generalized rates in first order in Tjy 

X2X2 

and ip(x) is the digamma function. It is, thus, the isospin 
generation term , describing combined Andreev pro- 
cesses that involve both the normal and the supercon- 
ducting leads, which reduce Q/(—e) below 1. 



V. CONCLUSIONS 

The theoretical description of sub-gap (Josephson and 
Andreev) transport through quantum dots with strong 
Coulomb interaction coupled to normal and supercon- 
ducting leads simplifies substantially in the limit of 
infinitely-large superconducting gap A and weak tunnel 
coupling Tn to the normal lead. For experimental de- 
vices with finite A and larger Tn these calculations may 
still be used as an approximation. The assessment of the 
quality of this approximation to describe sub-gap trans- 
port has been the focus of this paper. In particular we 
came to the following conclusions. The position of the 
ABSs strongly depend on A, i.e., a A — >• 00 approxima- 
tion is insufficient. These ABSs define threshold volt- 
ages in the current-voltage characteristics at which new 
transport channels open. The numerical values for the 
Josephson and Andreev current between these threshold 
voltages, however, seem to be nicely approximated by the 
A — » 00 limit. 

Similarly, we find that the next-to-leading order correc- 
tion in the tunnel coupling Tn to the normal lead mainly 
yields small quantitative corrections. But there are also 
qualitative differences: the position of the — tt transi- 
tion in the Josephson current and the peak positions in 
the Andreev current are shifted. 

In conclusion, we find that the A — >• 00 and weak- 
coupling calculations provide a very good approximation 
for sub-gap transport with the limitation that the posi- 
tions of the ABSs for finite A are not properly described. 
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To address the latter, we proposed a resummation ap- 
proach that substantially improves over mean-field treat- 
ments and favorably compares with more elaborate NRG 
calculations. 
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